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ABSTRACT 


An analytical model of a jet injected normally from a 
flat plate into a uniform crossing flow was modified to 
provide an improved method of predicting the interference 
effects arising from the complex flow fields induced by ship 
bow-thrusters. This model is an extension of previous work 
based upon a description of the jet as a series of distrib- 
uted vortices. The analysis takes into account the position 
of the effective source of the jet and the blockage due to 
the presence of the jet in the crossflow. Improvement in 
the predicted pressure distribution on the plate was obtain- 
ed by adopting a new formula for determining the effective 
Jet source which more closely approximates the experimental 
results. Further improvement resulted from the incorpora- 
tion of a symmetric foil shape to simulate the presence of 
the jet and its highly turbulent wake. The accuracy of the 
model was evaluated for representative jet-to-crossflow 
velocity ratios. Good agreement with experimental results 
was achieved for large portions of the interaction field 
with particular improvement in the near-field and down- 


stream regions. 
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I. INTRODUCTION 


Since man first set sail, ships and boats have been 
steered by stern-mounted rudders or rudder-like devices. 
The bow has been left relatively uncontrolled. As the size 
of ships increased, they became more difficult to maneuver 
in close quarters, as when entering port. Tug boats were 
often required to provide lateral thrust at both ends of the 
ship. This was a satisfactory solution to the problem for 
many years. It did not provide the degree of independence 
desired by ship operators, however, and other methods were 
considered. Although experimental lateral thrust units were 
used in the mid-nineteenth century, the first practical unit 
was not installed until 1955. This unit, a through-tunnel 
type bow thruster installed in the "Princess of Vancouver" 
proved highly successful [1]. Since then, bow Pheeeers 
have steadily gained in popularity. The advantages of these 
units include providing a greater degree of maneuverability, 
increased safety, increased ability to operate in rough 
weather at oil rigs and research sites, and reduced tug 
charges. 

In the past twenty-five years, ports have become increas- 
ingly crowded, labor and fuel costs have greatly increased 
tug boat charges, and ships have been built that are much 


larger and less maneuverable. Thus lateral thrust units 





have not only become economical, but they are standard 
equipment for many new~construction ships. 

As more bow thrusters came into use, it became apparent 
that although they work well at zero or very low forward 
speed, their effectiveness drops considerably as the ship's 
speed is increased [2]. In view of the large number of 
these units currently installed, there is Suprisingly little 
information available concerning thruster/hull interactions. 

Soon after the first unit was installed in 1955, it was 
learned that the operation of a tunnel-type thrust unit 
causes a low pressure area to form on the ship's hull [2] 
(see Fig. 1). To find the net resultant thrust, this suc- 
tion force must be integrated over the effected area and 
subtracted from the thrust of the unit. When a thrust unit 
is operated, water surrounding the jet is entrained, as in 
an eductor, cauSing a low pressure area to be formed on the 
Ship's hull in a small area about the jet. In addition, if 
a ship 1s moving forward, a plume is developed aft of the 
thrust unit and a system of vortices is formed. These vor- 
tices induce a velocity in the water at the ship's hull, 
causing a pressure drop. Additionally, the jet itself acts 
aS an obstruction to the flow of water past the ship, caus- 
ing the flow to change speed as it passes, thus causing 
another pressure varlation. These induced velocity and 


blockage effects are so severe that a thrust unit's 
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effectiveness can drop as much as fifty percent with a speed 


increase of only three knots [1]. 


The effects of bow thruster jet/hull interactions upon 


a ship motion can be better understood by considering the 


equations of motion for a ship moving in a horizontal plane. 


The general forms of the equations are [3]: 


m(v+u) 


m(u-vwy) = X Surge equation 


4 Sway equation 


>Y = N Yaw equation 


The linearized forms of these equations are sufficient 


for this analysis and, making the ship's center of gravity 


the coordinate system reference point, they become 


respectively: 


where: 


-Xfu-u,) + (m-x.)U = X (1) 
MT + (m-Y.)v = (Y-mu,) © = Y= =y (2) 
Sy - NV -- Ne + (il AN a) se = N (3) 


mass of the ship 

velocity in the x-direction (initial value us) 
velocity in the y-direction 

yaw angle 


angular velocity = yp 


mass moment of inertia 


Je 





Ko-sesxeitaeton forces ian the x-direction 
Y = excitation forces in the y-direction 
N = excitation torque 


Dots indicate differentiation with respect to time. 


The subscripts u, 0, v, v, xr, xr, denote differenti- 
ation with respect to the subScript variable. 


The left sides of these equations contain the velocity 
and acceleration-dependent terms such as added mass, added 
inertia, and hydrodynamic drag. The right sides implicitly 
represent the total forces in the x and y directions and the 
total moment about a vertical axis through the center of 
gravity. Although these terms are used to represent various 
excitation forces acting on a ship due to rudders, stabi- 
lizers, propellors, etc., this analysis will deal solely 
with the forces and moments caused by a through-tunnel 
bow=mounted thrust unit. 

In the y-direction, the force due to a bow thruster may 


be written as: 


iE S 
where: 
Y, = “Mus, the ideal jet thrust 
Y, = the induced isuctmienmzerneceonethe hull 
us = the mass flow rate of the jet 
us = the jet velocity (positive to starboard). 


ILA 





Inspection of Fig. 2 reveals that “Y/Y, 1s nearly linear 
in the pu / a. | Bueko range Of 21 tO .3. Using this velocity 
ratio range aS an example, an approximation of the right 
Side of equations (2) and (3) can be made. This portion of 
the curves is of interest because the largest decrease in 
thruster performance takes place here. Since the slope and 
intercept of this line are dependent on hull shape, an 
average performance curve has been plotted in Fig. 3 to 
represent a "typical" hull. A straight line having a slope 
of 5.3 and an intercept of -0.65 is used to approximate the 
curve in the ju /a, | BAbTOu(Menbeatiter R) of .1- to .3. 

Thus: 


pata —oeok = 60 


Then since 


Vie =... 
t i 
Ve Met (oR =. 65 4 
Z gia ) (4) 
and 
Y= greancl. 65 -5.3R) ile re Se (5) 


By substituting this equation into equation (2), the 
motion in the y-direction is explicity coupled to the motion 
in the x-direction through the dependency upon R. 

The total turning moment N is the moment produced by 


the force Y acting at a distance ve from the center of gravity. 


ae 





This characteristic moment arm varies non-linearly as the 


ship's speed changes and is given as: 


a XY, te on oe 
c ne 
where: 
eo the distance from the center of gravity 
(ey eleva) we gue cm nbbe lite 
aoe the distance from the center of gravity 


to the center of action of the hull 


Suction force. 


Chislett and Bjorheden [{2] found that < varies linearly 
with ship's speed. The curve in Fig. 4, which is adapted 
from their data, has a slope of -.7 and an intercept of 1.0. 


Therefore: 


~*~ 
iT 
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i 
~J 
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NO 


(6) 


When equations (4) and (5) are substituted for Yor eran Claw Ye, 


t 


the location of the center of action of the combined force 


system is given by: 


_ (3.71R* - 5.76R + 1.65), 


Gl Gs) 5 Glee ie 


and 
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Thus the yaw equation is coupled to the sway and surge 
equations. Note that equation (7) is only valid for velocity 
ratios between .1 and .3 for this particular fictitious hull. 

Although this analysis has been greatly simplified, it 
demonstrates the complication of the relationships between 
turning moments and ship's speed when a lateral thrust unit 
is used. It can be seen from Fig. 2 that model tests of 
each ship type are necesSary in order to accurately predict 
the performance characteristics of a lateral thrust unit. 
With a thorough understanding of the fluid mechanics of bow 
thruster systems, thrust units could be designed to reduce 
the hull suction effect and thus make them more usable under 
a wider range of conditions. An accurate analytical method 
for predicting jet/hull interactions would be useful in 
supplying this information without lengthy and costly model 
tests. Unfortunately, an accurate model that can be super- 
posed with other effects is not available. It is the purpose 
of this thesis to propose an analytical method for predicting 
the jet/hull interactions of a through-tunnel tyve lateral 


Garust unit. 
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II. BACKGROUND 


As demonstrated in the preceeding section, the ship's 
hull/thruster interactions, when a ship has way on, can 
seriously degrade a bow<-thruster's effectiveness. In order 
for naval architects to design a bow thruster installation 
so that it can be used effectively when the ship 1S moving, 
a complete knowledge of these effects is required. Unfortu- 
nately, only a few waterborne tests have been performed. 
These experiments (Fig. 2), while clearly demonstrating the 
@merall effects on turning ability, do not provide informa- 
tion on the actual pressure distribution on the hull. For 


an incompressible flow, the jet path and the induced flow 
z 


Lo @ 


are dependent mainly on the momentum flux ratio p40; /0,U 


and are relatively independent of the Reynolds number [4]. 
Thus an air jet injected into a crossing air flow will induce 
a reaction similar to that of the water jet of a thrust unit 
at the same velocity ratio. There have been a number of 
studies of air jets emitted from a flat plate into a crossing 
air flow in connection with the design of VSTOL (Vertical and 
Short Take Off and Landing) aircraft. Notable studies have 
been made by Bradbury and Woods [4] and Fearn and Weston [5]. 

While the results of these experimental studies are use- 
ful to the naval architect when making an initial design 


intended to improve thruster efficiency, they provide no 
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method for evaluating the resultant configuration. Lengthy 
and expensive model or shipboard tests are required for this 
purpose. An accurate and versatile analytical model of the 
flow interaction effects would provide a quick and inexpen- 
Sive means of evaluating new designs, thereby reducing the 
required number of model tests. Some empirical and semi- 
empirical models have been developed. While they can predict 
the pressure distributions on a plate, they cannot be easily 
Superposed with other effects such as a secondary jet. 
Because a fluid jet injected into a crossing stream 1S a 
turbulent, three dimensional, highly non-linear problem, 
empiricism to some degree has been used in all models to 
date. In this study an effort has been made to eliminate or 
at least reduce the empirical inputs to existing analytical 
models while improving their accuracy. 

A review of available models revealed that Wooler's formu- 
lation [6] minimized empiricism, requiring only an experi- 
mentally determined jet trajectory. Hence his model, as 
modified by Waterman [7], was selected as a basis for further 
investigation. 

In order to visualize the interaction field and determine 
the accuracy of an improved model, once formulated, accurate 
experimental results were required. A review of available 


studies revealed that the pressure contours* for a given 


* 
In thisthesis, values of local pressures, P, are expres- 
sed in the form of pressure coefficients, Cor where: 
Pp-P 


oe Ly 





velocity ratio vary considerably from one report to the next, 
as illustrated in Fig. 5 and 6. This difference is especially 
noticeable in the downstream region for values of }U5/U,, | 
(hereafter VR) of about 4, and 1S present to some degree at 
all velocity ratios. Figure 5 seems to imply that the induced 
flow may not be solely a function of velocity ratio since the 
different studies were not made at the same freestream 
velocities. However, considering studies made with the cross- 
flow velocity nearly equal, aS in Fig. 6, it is seen that the 
discrevancies are more probably due to variations in the 
equipment and test procedures used. This serves to point out 
the uncertainty of the available experimental results. The 
results of Fearn and Weston [5] agree well with those of 
other researchers over a broad range of velocity ratios. 
Additionally their report was the most complete of those 
reviewed. It contains not only the usual plots of pressure 
contours, but also raw data, time variation of pressure 
information, and a comparison with other reports. For these 
reasons, the results reported in [5] were selected as the 
basis for determining the accuracy of the analytical model. 
Before attempting to develop an improved model, a thor- 
ough review of Wooler's and Waterman's work was made in 
order to understand their aporoaches and verify their results. 
A rederivation of these basic formulations will not be pre- 
sented here as it can be found in references [6] and [7]. It 


was noted that Waterman's adoption of a conforming vortex 
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model was a considerable improvement to Wooler's formulation. 
In the conforming vortex model the vortices are constrained 
to lie within the jet rather than being tangential to the jet 
axis as in Wooler'ts model. By making this alteration, the 
horseshoe vortices are made to conform to the curvature of 
the jet described by the empirical trajectory equation (see 


fey. 7): 


= B[cosh(Z/Bd) - 1} (8) 


ul x 


Where: 
x = the distance down-stream from the jet center. 
d = the jet diameter. 
z = height of the jet axis above the plate. 


B = trajectory parameter = 19 (VR)°. 


The interference velocity at an arbitrary point on the 
plate is determined by applying the Biot-Savart Law. The 
trailing vortices are approximated as straight line segments. 
As these segments become small, the linear approximation can 
. be made to conform to the curvature of the jet. The total 
interference component at a particular point due to the 
entire jet results from the summation of all the elemental 
@entributions. 

A comparison of the results of both formulations reveals 
that Wooler's work more closely agrees with the experimental 


findings of Bradbury and Woods [4] than does Waterman's, 


i 





even though the trailing vortices in Wooler's model do not 
follow the jet trajectory as in the actual case. This result 
is, in part, due to the fact that in Wooler's formulation the 
perturbation velocity at a particular point on the plate is 
caused by vortices acting at different distances from the 
plate. In Waterman's model the interference velocities 
result from the trailing horseshoe vortex contributions act- 
ing at a Single radius. Thus the interference velocities 
predicted by Waterman's conforming vortex model are generally 
higher and hence produce lower coefficients of pressure than 
those of Wooler. In addition to modifying Wooler's model, 
and noting studies by Keffer and Baines [8] and Schmidt [9], 
Waterman introduced the concept of an effective origin of the 
vortex sheet at a distance, Ze above the plate given by 


the furmula: 


_ VRVVR 
eee eS) eh, 


Applying this expression to the conforming vortex model, 
Waterman produced the same, if not better far-field agreement 
(Fig. 8). However, the model lacked accuracy in the near- 
field. 

Waterman further noted that the discrepancies observed in 
the upstream region can be at least partially attributed to 
the blockage of the freestream flow by the jet. It was noted 


that the downstream area is dominated by separation, vortex 
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shedding, and turbulent entrainment. Waterman interpreted 
the blockage effect as that of a circular cylinder in an 
irrotational flow. Using potential flow theory to predict 
the velocity perturbations in the uniform flow caused by a 
doublet centered at the jet origin, Waterman, after some 
numerical experimentation, arrived at an estimate for an 
effective jet diameter given by .96V¥YVR . The results of 
this model (Fig. 9) are significantly better than Wooler's 
model in the upstream region and far field. This improve- 
ment was accomplished at the expense of introducing two 
additional empirical relationships, the position of the 


effective source and the effective jet diameter. 


Zn. 





ct. Sie ReyVEMENT IN THE EXISTING MODEL 


Evaluation of Waterman's "fully-corrected conforming vor- 
tex model" (Figs. 9 and 10) revealed that the correlation of 
the analytical model with the experimental data in the area 
just forward of the jet is adequate at each of the velocity 
ratios reported. However, his model is accurate in the far- 
field only for the velocity ratio of 8.0 and in the near- 
field to the sides of the jet for a velocity ratio of 4.0. 
The correlation is poor in the area aft of the jet for all 
velocity ratios. It was apparent that the near and mid- 
field regions, especially aft of the jet, required the 
greatest attention. 

In order to gain a better understanding of the effects of 
the empirical inputs to Waterman's model, it was run a num- 
ber of times uSing various values of the constants in the 
trajectory, effective source, and effective diameter 
equations. It was learned that small changes in the trajec- 
tory parameter, B in equation (8), cause little change in 
the overall pressure distribution. The correlation with 
experimental data was improved at some velocity ratios while 
it diverged at others for a particular value of B. A value 
of B = .19(VR)~, as before, resulted in the best overall 
correlation of the pressure distribution on the plate. 


Changing the effective diameter of the jet causes 
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major changes in the near-field pressure distribution, but 
has little effect in the far-field. This was expected con- 
Sidering the short range of the disturbance caused by a 
doublet in a uniform flow using inviscid relationships. The 
model proved to be particularly sensitive in both the near 
and far-field regions to changes in the position of the 
effective source. Since the quality of the model depends 

SO greatly on the effective source correction, the equation 
for it was investigated more closely. 

Equation (9), the empirical formula for the effective 
source height, was formulated by Schmidt [9] who refers to 
Jordinson [10] and Keffer and Baines [8] for the experimental 
data used as a basis for his work. The Jordinson paper could 
not be obtained for review. The data presented by Keffer 
and Baines, however, does not fit the equation presented by 
Bemmiadt. The equation, Z/4 = .8VVR fits the Keffer and 


Baines data much better as seen in Table Tf. 


TABLE I 


COMPARISON OF VALUES OF Z_/d 


Experimental 
Keffer & 
VR Schmidt Baines .8V¥ VR 
4 53 Peo 160 
6 i 1.95 1.93 
8 digpadovl li PAP EY DOS 
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However, Z,/d can vary from .54 We Reco .¢ % VR and still be 
in the range of the heights given in Table I. Using an 
effective source location of Z_,/4 = .8VVR in the Waterman 
model significantly improves the correlation. But, if 
py = .6VVR is used, the correlation is much better even 
without the effective diameter correction. Thus, by estab- 
lishing the effective source at .6 VVR diameters above the 
plate, the full jet diameter can be used in the model and 
the diameter correction is eliminated. See Fig. 11 for an 
example of the improved correlation resulting from this 
change. 

Concentrating on the area aft of the jet, 1t can be seen 
PP ecomparing’ Fig. 8 with Fig. 10 that this area is primarily 
affected by the blockage near the plate due to the jet. In 
making multiple runs of Waterman's model and varying the 
parameters of the vortices within the plume, it was confirmed 
that they have little effect on the near-field pressures. 
The far-field to either side of the jet is dominated by the 
vortices within the plume. Therefore, to improve the model's 
- accuracy in the near-field, the analysis centered on the jet 
blockage effect. Wu and Wright [11] proposed that a 
Fictitious afterbody exists downstream of the jet. fThis 
afterbody encompasses the highly turbulent, three-dimensional 
flow caused by vortex shedding and entrainment of fluid 
within and downstream of the jet. Outside of this area, the 
fluid streamlines are fairly smooth and close to those of an 


inviscid flow. With this concept in mind, consideration was 
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given to using an elongated shape such as an ellipse or a 
symmetrical foil to approximate the blockage effect. 

In order to establish a probable shape for the afterbody, 
the coefficients of pressure along the corssflow axis from 
the forward edge of the jet to seven diameters aft of the 
jet were plotted using the data of [5] (Fig. 12). These 
pressure contours resemble those in reference [12] and [13] 
for a symmetrical foil in a uniform flow. Thus a symmetrical 
foil was selected to approximate the jet and the afterbody 
near the plate. 

To determine a length for the foil, it was assumed that 
the foil should end where the coefficient of vressure reaches 
an undisturbed (zero) value. It was not expected that the 
coefficient of pressure should necessarily rise above zero 
aft of the jet due to severe dissipation in the near-wake 
region and the interference effects of the vortices within 
Geemorume. Indeed the coefficient of pressure leveled-off 
at approximately -0.1 for all available data. Therefore, 
-0.1 was selected as a first approximation of the value of 
the coefficient of pressure used to determine the length of 
the foil. The corresponding distances from the jet center 
to the trailing edge were determined to be 6,2, 3.0, and 2.5 
fetceciameters for the velocity ratios of 3.9, 8.0, and 10.0, 
respectively. 

For simplicity, a standard Joukowsky foil was selected 


Since it can be easily approximated through mapping methods. 
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Referring to Fig. 13, the Joukowsky transformation, 


is used to transform the foil in the z-plane to a cylinder 
in the ¢-plane so that the foil can be treated as a doublet 
in a uniform stream. The perturbation velocities due to the 


presence of the foil are then: 


and 


Where: 


u = the surface perturbation velocity in the direction 
of the cross-flow due to the presence of the jet 


and afterbody as modeled by the foil. 


v = the surface perturbation velocity pervendicular to 
the cross-flow due to the presence of the jet and 


afterbody as modeled by the foil. 


(=) Se ene distance from the forward focus of the 
foil (corresponding to the point (-%,0) in the 


C-plane) to the trailing edge. 


a = the radius of the transformed circle. This term 
1s a function of r, the distance from the forward 
focus of the foil to the leading edge. 
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Appendix A contains a complete derivation of the above 
transformation. 

The thickness of the foil is a function of both r and 2. 
As described above, the foil length is a function of the 


velocity ratio so that the length & is given by: 


ie! CEO} 


Initially it was assumed that the forward focus of the 
foil should be at the center of the jet. This makes r equal 
to one-half the jet diameter. When substituted for the 
doublet in the fully corrected model, this method provided 
Significant improvement in the downstream area. The up- 
Stream correlation was poor, however, since the foil was too 
thick and too blunt at the leading edge. It was necessary 
to adjust the value of r to provide a thinner and more 
Streamlined body. Next a foil with its thickness equal to 


one jet diameter was tested. Using the formula, maximum 


thickness/d = 5.2(m/2 - Pima [13], where m is a function 
of r. It was found that if r/d = .12d/2, where 2 is given by 
equation (10), the thickness is nearly equal to the jet 


diameter for the entire range of velocity ratios. The foil 
is then moved forward so that its leading edge coincides 
with that of the jet. This configuration (App. B) produced 
the results shown in Figs. 14, 15, and 16 in comparison with 


experimental values [5] for velocity ratios of 4.0, 8.0, and 
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10.0 respectively. Radial and circumferential pressure 
distributions for velocity ratios of 3.9 and 8.0 are plotted 
against the experimental findings of Fearn and Weston [5] 
Baerigs. C-l1 through C-52. 

Attempts were made to use other fictitious shapes such 
as ellipses and ovals as well as distributed source-sink 
combinations. The foil shape, however, produced the best 
results. So the simple relationships for the near-field 
blockage effect given above were retained in the final 
model. 

Reviewing the results, it is seen that for a velocity 
ratio of 8.0, the correlation between this model and the 
experimental data is very close except in the area of 
highly turbulent flow immediately aft (within two diameters) 
of the jet. This is a region of uncertainty for all velocity 
ratios due to the three-dimensional flow caused by entrain- 
ment, vortex shedding, and other highly-viscous wake 
phenomena. Since this area is small, it contributes little 
to the overall thruster/hull interaction. For a velocity 
ratio of 3.9, the correlation is excellent except for the 
mid- and far-field regions aft of the jet centerline. 

This was expected since it was determined earlier that 
the blockage effect of the jet has little effect in the 
far-field. It was observed that if r is such that it causes 
the foil to have a circular shape at the leading edge the 


correlation in the area forward of the jet is improved. 


28 





However, the agreement deteriorated in the areas to the 

Sides and aft of the jet. By varying r/d as a function of 
the location of the point at which the velocity perturbations 
are being calculated, the shape of the foil can be altered 
so that it is nearly circular at the leading edge and becomes 
thin aft of the jet centerline. This configuration worked 
better, but it failed to significantly improve the overall 
pressure distribution resulting from r/d = .12d/2%. This 
observation did lead to the idea of uSing a half-cylinder 
Shape for the jet blockage. That is, the fluid forward of 
jet "sees" a cylinder but once past the jet centerline 

the influence of the cylinder is rapidly diminished. The 
results of this configuration are presented in Fig. 17 from 
which it can be seen that while this is a Significant 
improvement over Waterman's model, it is not as accurate 

as the full foill-shape model. It is concluded that further 
improvement of the model must involve a more sophisticated 
treatment of the vortices within the plume as well as the 


blockage effect of the jet. 
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IV. CONCLUSIONS 


Cemecvrdencea by Bags. 24, 15, 16, and C-1 through C~52 
this analysis adequately predicts the pressure distribution 
on a flat plate from which a jet is emitted into a crossflow. 
Application of the model is limited, however, to velocity 
ratios VR = 3.9 or greater. This is due to the method by 
which the length of the fictitious afterbody 1s determined. 
Above velocity ratios of about 4, this length decreases for 
increasing VR, but plots of the coefficients of pressure 
along the crossflow axis Similar to Fig. 12 indicate that 
the afterbody length should decrease also for values of VR 
less than about 4. This is consistent with the observations 
of Fearn and Weston. They reported that an anomaly in the 
pressure distribution in the wake occurs at a velocity ratio 
oud 3.9. The implication is that the rate of pressure recovery 
in the wake region is a minimum for a velocity ratio of about 
nee, (>). 

No attempt has been made to accurately predict the pres- 
Sure distribution in the highly turbulent wake of the jet. 
This area, having a width of about one jet diameter, extends 
downstream from the jet along the crossflow axis. The pre- 
dicted values of the plate pressures reported here cannot be 
relied upon since most of this area is contained within the 


fictitious afterbody. 
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Figure 18 shows the results of this analysis plotted as 
the ratio of the induced side force to the jet thrust, 
me e/ Ye versus the velocity ratio VR. These results are 
compared with NPL pressure measurements on a flat plate and 
the model tests of English [1] and Chislett and Bjorheden [2]. 
The slope of the predicted ~¥6/Y, curve is approximately 5.5, 
whereas that of the experimental data is about 5.3 for both 
the flat plate and the model tests. While this analysis 
predicts somewhat larger values for the overall (integrated) 
interference effects, the small difference between these 
slopes indicates that the nature of the relationship between 
the interference effects and the velocity ratio predicted by 
the analytical model is correct. Fearn and Weston did not 
report overall interference effects so the correlation between 
their work and the NPL experiments is not known. Since the 
assumptions concerning the shape and size of the afterbody 
were based on Fearn and Weston's data, the degree to which 
this analysis is in error cannot be determined from Fig. 18 
alone. Even with these deficiencies, it is felt that this 
analytical model adequately predicts the pressure distribu- 
tion on a flat plate and can be used as a first approximation 


for design purposes. 
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Vo PRECemiIrNDalLONsS ee OR FURTHER STUDY 


In Chapter III it was noted that afterbody shapes other 
than a symmetric foil have been tried without significant 
improvement. This does not mean that this approach showed 
no promise for possible improvement in the model, but only 
that the actual shapes tested were not the best. It is felt 
that this model can be improved by changing the shape used 
for the jet and afterbody. 

As reported earlier, the crossflow fluid particles near 
the plate apparently see a cylindrical shape when approaching 
the jet. As they pass the jet, some particles are entrained, 
some join the turbulent flow within the wake, and, further 
out, they pass as if the jet and wake together are a solid 
object. Just downstream of the jet, along the crossflow 
axis, there is an area of high entrainment. This causes a 
high local pressure drop out to about one jet diameter for 
VR = 8. Pressure recovery 1S rapid past this point. This 


leads to the idea of uSing a shape such as: 


This shape has a cylindrical leading edge, the maximum 


thickness is just downstream of the jet center, and then it 


ee 





mnarrows quickly to a relatively thin tail. The total length 
is the same as that for the symmetric foil. 

It is expected that this shape can be formed using the 
Joukowsky transform of the existing model by varying r as a 
Bamction of the angular position of the particular point at 
which the perturbation velocities are being calculated. 
Although a functional relationship which will work over a 
range of velocity ratios has not yet been found, it 1S pos- 
sible that one exists. This idea might also be pursued 
using panel (distributed singularity) methods. 

Another area of uncertainty during the course of this 
study was the empirical jet trajectory equation. Reports of 
experiments to determine the jet axis provide a wide range 
of possible trajectories. If the trajectory equation used 
in this model is accepted as correct, the trajectory con- 
stant B = 19VR? still may not be accurate over a wide range 
@eevelocity ratios. In particular, it is suspected that it 
may be in error for velocity ratios near 4 where an anomaly 
in the plate pressure distributions occurs. The cause of 
this anomaly is not clear, but it is suspected that this is 
the point where the jet momentum is so low compared to that 
of the crossflow that the jet reattaches to the olate within 
a few jet diameters. If so, the jet trajectory would be 
quite different than that which is currently modeled. 

Finally, the entrainment of the crossflow fluid into the 


jet is implicity accounted-for in the formulation of the 


25 





horseshoe vortices. However, the pressure drop associated 
with the jet entraining the surrounding fluid is not explic- 
itly determined. This model predicts a zero pressure distri- 
bution at very high velocity ratios (zero ship speed). 
Further improvement of this model in accounting for entrain- 
ment is needed if it is to be used to predict the performance 


of bow-thrusters in ships with no way on. 
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Figure 7. Jet axis and conforming vortex loop system. [7] 
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APPENDIX A 


DERIVATION OF THE JOUKOWSKY TRANSFORMATION 


The Joukowsky transformation, 2 = ¢ + 22 /z, 1s used to 
transform a foil in the z-plane to a cylinder in the ¢-plane 
so that the foil can be treated as a doublet in a uniform 


flow. The terms used in this transformation are: 


Z=- the complex position vector of a point in the 
Z-plane Z = x + ly 

c - The complex position vector of a point in the 
t-plane. ¢ = & + in 

2 —- A positive real constant equal to 1/4 the distance 
between the foci of the foil. 

a - The diameter of the cylinder in the c'-plane. 

m - The distance between the yn and 7» axes. This is a 
BUNS Lon Of ir. 

r - The distance from the forward focus to the leading 
edge of the foil. That is r = |/c - 22] where -c is 
the abscissa of the leading edge. 

c'-plane - The linearly transformed c-plane used to give 


Pe miorleentchitess, 6 = oe ey 


From the Joukowsky transformation and Fig. 13 it is seen 
einat : 


ae 252 St Z~ - 42° ) (Al) 
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The complex potential of a cylinder in a uniform stream 


1s given by: 


mee) Slee tl a/c: ) (A4) 


The complex velocity is then found as follows: 


=o = dw = dw dc! ug 
(8 Se ul Vy az dc! dt az, 
Gina ) dice az 
az > ae ge’ / ee Ci 
then substituting 
2 12 2) 
du Ly (ie a) 
az, co : 
G es 
Bema Since c -a=c«-2 
2 t 
dz © (= eS) (A6) 
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